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1.1  Basic  Assumption 

We  assume  that  the  shell  i3  a rigid  body.  The  configuration 
at  any  tine  may  therefore  be  defined  by  the  rectangular  coordinates 
of  the  center  of  mass  and  three  Eulerian  angles  giving  the  orienta- 
tion of  a set  of  axes  moving  with  the  shell.  The  state  of  motion 
of  the  shell  may  be  defined  by  the  time  rates  of  change  of  these 
six  generalized  coordinates,  making  together  with  them  tvielve 
variables  capable  of  independent  variation. 

We  assume  that  the  system  of  forces  and  torques  acting  on 
the  shell,  consisting  of  its  weight,  and  of  aerodynamic  forces 
and  torques,  depends  only  on  thsse  twelve  variables,  except, 
perhaps,  for  small  correction  terms;  more  detailed  assumptions 
will  be  made  later. 

We  assume  that  the  shell  has  an  axis  of  dynamical  symmetry 
passing  through  its  center  of  mass. 

For  simplicity  in  the  following  discussion  we  also  neglect 
the  rotation  and  curvature  of  the  earth,  and  suppose  that  there 
is  no  wind. 


1.2  The  Variables  Used  to  Specify  the  Motion 

We  use  twelve  variables,  nine  of  which  differ  from  the 
above  variables. 

Firstly  we  take  three  rectangular  coordinates  x,  y,  and  z, 
of  the  center  of  mass  of  the  shell  in  a fixed  coordinate  system; 
the  x-direction  horizontally  forwara,  the  y-direction  vertically 
up.  and  the  z-direction  horizontally  to  the  right.  The  plane 
x 0 y contains  the  line  of  fire. 
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v,a,6  Secondly,  we  take  the  speed  v of  the  center  of  mass  and  the 

azimuth  a,  measured  from  the  x-direction  to  the  loft,  and.  angle 
of  elevation  0,  measured  from  the  horizontal!,  upwards,  of  the  di- 
01  rection  01  of  the  velocity  of  the  center  of  mass,  so  that  we  have 

x * v cos  0 cos  a ^ 
y * v sin  0 L 

z = -v  cos  0 sin  a 

to  relate  the  tine  rates  of  change  of  the  coordinates  x,  y,  and 
z,  to  these  new  variables. 

(0  is  not  any  particular  point,  but  is  an  origin  of  reference  for 
directions,  which  will  be  represented  diagrammatically  by 
points  on  a sphere  of  center  0.) 

ON  ON  is  to  be  the  direction  at  angle  a forward  from  Oz,  in  a 

horizontal  plane,  and  about  which  0 is  measured. 


Next  complete  the  determination  of  the  rectangular  set  of 
02,03  directions  01,  02,  and  03,  by  making  the  plane  102  contain  the 
OA  direction  OA  of  the  axis  of  dynamical  symmetry  of  the  shell. 

Then  102  is  the  plane  of  yaw,  and  the  angle  0 from  the  vertical 
plane  yOl  to  the  plane  201,  is  0 the  angle  of  rotation  of  the 
plane  of  yaw,  measured  in  the  righthand  direction  about  01. 

01,  02,  and  03,  is  then  a set  of  rectangular  axes  with  orienta 
tion  in  terms  of  the  fixed  directions  Ox,  Oy,  and  Oz,  given  by 
0 angles  a,  0,.and  0.  Resolving  the  corresponding  angular  velocities 

a about  Oy,  0 about  ON,  and  0 about  01,  in  the  directions  01,  02, 
and  03,  we  obtain  the  angular  velocity  components  and 

of  this  frame  of  reference  referred  to  itself,  namely, 

^ + a sin  0 

$2  * i cos  0 cos  0 + 0 sin  0 

■ -a  cos  9 sin  0+0  cos  0 
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We  shall  use  the  letter  A for  the  principal  moment  of 
inertia  of  the  shell  about  the  axis  of  dynamical  symmetry. 

On  account  of  the  symmetry,  the  moments  o f inertia  about  axes 
perpendicular  to  OA  through  the  center  of  mass  are  the  same; 
they  will  be  denoted  by  B. 

6,  the  angle  of  yaw  is  the  angle  10A;  we  shall  call  the 
direction  perpendicular  to  OA  in  the  plane  10 A,  (which  includes 
02),  OB.  The  principal  moments  of  inertia  of  the  shell  about 
the  directions  OA,  OB,  and  03,  are  then  A,  B,  and  B. 

ifr  is  to  be  the  angle  of  rotation  from  the  plane  AOB  to  a 
plane  moving  with  the  shell. 

We  take  the  angles  0,  6,  and  ifT,  as  the  third  3et  of  vari- 
ables to  specify  the  motion.  They  are  Eulerian  angles  giving 
the  orientation  of  the  shell  relative  to  the  tangent  direction 
01  to  the  trajectory  of  the  center  of  mass  and  the  vertical 
plane  lOy  through  01. 


B 


2 


$>A  ,f 
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For  the  fourth  set  of  three  variables  we  take  3^/^and 
the  components  of  angular  momentum  of  the  shell  about  01,  03, 
and  OA.  The  components  of  angular  velocity  of  the  shell  about 
OB,  03,  and  OA,  are  obtained  by  re  solving  $■-  and  in 

these  directions,  and  taking  account  of  the  angular  velocities 
corresponding  to  6 and  ijr,  giving 

■$2  cos  6 - ^ sin  6 

^3+5 

A sin 

The  corresponding  components  of  angular  momentum, 22,  say, 

A , and  ijr,  are 

22  = B ( $2  cos  6 - 2/^  sin  6) 

A=  B (zPjj  + 6) 

A (zPg  sin  6+2^  cos  6 +j^) 
so  that  £ , the  component  about  01  is  given  by 

£=  A cos  6(^2  sin  6 + ^ cos  6 +jfr)  - B sin  6($>  cos  6-2^  sin  6). 

The  variables  x,  y,  z,  v,  c,  0,  0,  6,^",  <£>  , A , and  jjp',  are 
all  given  if  the  configuration  and  velocity  of  the  shell  are 
given.  This  would  not  be  true  if  we  tried  to  use  0,  6,  and 
y',  for  these  would  involve  a and  0,  and  therefore  x,  y,  and 


The  component  of  angular  momentum  ir.  the  direction  02  is 
B cos  6(^2  cos  6 “2^  sin  6)  + A sin  6(2P-2sin  6 +2^  cos  6 + ^) 

ijr-  £ cos  5 f 
e sin  6 

j7*  cos  6 ~ $ 

while  22*  sin  6 
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1.3  The  Equations  of  Motion 

The  absolute  rate  of  change  of  a vector  whose  components  in 

the  01,  02,  and  03  directions  are  L,  L„  and  L-  has  components 

■*•**;>  > 

H * ^2L3"  ^3  L2 


and 


M 


If  M is  the  mass,  the  linear  momentum  has  components  in 
these  directions  Mv,  0,  and  we  have  the  equations 


M v = F, 


V*.FJ 


2^3  H v * Fj 
- M v •=  F3 


for  the  rates  of  change  of  momentum  equated  to  the  components  of 
force  acting  on  the  shell  in  the  01,  02,  and  03  directions. 

In  the  same  way,  the  conponents  of  angular  momentum  being 
3 ~ and  A , we  obtain 


sin  6 


C03  6 \ 


G G 

1,  2, 


'2  “ $ 003  ^ \ - G 

sin  6 I 1 

cRT\  sin  6 j 3?  . 1 2 

K*  (#~i^6C0S  S ^ “ ^2  ^ “ °3 


in  terms  of  the  components  of  torque  about  its  center  of  mass 
on  the  shell  in  the  01,  02,  and  03  directions. 

It  is  usual,  however,  to  regard  the  torque  components  about 

the  directions  OA,  OB,  and  03,  namely  G.  Gq  aw*  G,  as 

o,  o, 

primary,  and  to  express  the  equations  of  motion  in  terns  of  them. 
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(*A  cos  6 +°2  sin  6 

GB  = "G1  sin  ® + G2  cos  6 
Combining  the  first  too  of  these  equations  vie  obtain 

sin  6 + G^  cos  6 * G^  , 


the  remaining  terras  cancelling,  which  equation  can  be  obtained 
directly  from  the  absolute  rate  of  change  of  angular  momentum 
about  OA. 


Eliminating  z9h  and  A we  find  from  the 


se 


equations  and  from  the  equations  of  the  last  section,  equations 
for  the  rates  of  change  of  our  twelve  variables  in  terms  of 
these  variables  and  of  the  components  of  force  F,  F0  and 

A, 

F^  and  of  torque  G^  and  G^. 


x « v coa  0 cos 
_y  » v sin  0 


-v  cos  8 sin  a 

F1 

rr~ 


e 


F F 

2 sin  0 3 cos  0 

“ W cos“0  ~Ev~  cosT 

F F 

2 cos  0-3  sin  0 

Hv~ 

F 

1 A 2 
B & W 

"if  cos  6 + F2  ^ y + *3  (cos  0 tan  q - cot 


r 

A 


B sin*1  6 

C03  6 


wr 


wr 


B sin  6 


( § - -jf  cos  5)  + u- 


F„ 


v sin  6 


- ( # - # cos  6)  ( 4 COS  6)  + G » (#  - cos  6)  F3 


B sin^  6 


‘ 2 “ Hv 
sin  6 
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Part  II.  The  Forces  and  Torques 


p?  f a 


2.1  The  force  and  Torque  Components 


The  force  system  acting  on  the  shell  is  made  up  of  gravi- 
tational forces  and  of  aerodynamic  forces  and  torques.  Gravi- 
tation contributes  terms  -M  g sin  $ , -M  g cos  & cos  0>  and 
M g cos  0 sin  0 to  F^,  Fg,  and  F^  respectively. 

L 


The  aerodynamic  forces  are  defined  relative  to  the  direc- 
tions 01,  02,  and  03,  but  are  taken  positive  in  their  expected 
directions,  the  drag  D in  the  opposite  direction  to  01,  the 
lift  or  cros3-wind  force  L in  the  direction  02,  and  the  Kagnus 
force  K in  the  opposite  direction  to  03.  For  a symmetrical 
projectile  rotating  from  OB  to  03  about  01  we  should  expect  a 
force  in  this  directionj  a ball  swerving  on  account  of  spin 
follows  its  nose. 

F,  =»  ~D  - mg  3 in  6 
F2  * L - mg  cos  0 cos  0 f 
F^  * -K  + mg  cos  0sin  0 
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Thus 


The  torques  are  defined  relative  to  the  directions  OA, 
OB,  and  03,  but  are  taken  positive  in  their  expected  direc- 
tions; the  spin  destroying  or  rolling  moment  I from  03  to  OB; 
the  yawing  or  overturning  moment  M from  OA  to  OB;  and  the 
Magnus  moment  J from  03  to  OA.  For  a projectile  with  Magnus 
force  in  the  opposite  direction  to  03  acting  in  front  of  the 
center  of  mass,  we  should  expect  a torque  in  this  direction. 

Thus  G.  **  - I 

A 

gb  - J r 

G3  « M 


2,2  Several  Assumptions  About  the  Aerodynamic  Forces  and  Torques 

We  assume  that  the  aerodynamic  forces  and  torques  depend 
only  on  the  properties  of  air  through  which  the  shell  is  pass- 
ing, which  we  shall  suppose  depends  only  on  y,  and  on  the  motion 
of  the  shell  relative  to  the  air;  that  is,  they  do  not  depend 
on  x,  2,  a,  9,  or  0.  Thus  the  aerodynamic  forces  and  torques 

, iS"  2?  A 

^depend  only  on  y,  v,  6,  jT , A p IT*  a manner  de- 

termined  by  the  properties  of  the  air,  specified  by  its  density 
p,  sound  velocity  c,  and  viscosity  n;  the  size  of  the  shell, 
specified  by  its  caliber  d,  and  its  dhape  relative  to  its  center 
of  mass. 


Here  we  have  tacitly  assumed  that  any  effects  of  lag  in 
the  distribution  of  airflow  can  be  taken  account  of  by  the 

dependence  on  f # , and  . 


Dynamical  similarity  then  requires  that  force  components 

2 2 2 3 

be  proportional  to  p v d and  torque  components  to  p v d , 

with  coefficients  depending  on  the  zero-dimensional  quantities 

v/c  =7 VI  s the  Mach  number,  • • * R,  'toe  Reynolds  number,  the 

number  ' p 

d.  ip  (j  ^ 

of  radians  turned  per  caliber  advanced,  — — , - ~ , and 

* v A ’ v B * 

~ , and  the  angles  6 and  . 


2.3  Aerodynamic  Forces  and  Tor:iues  for  a Shell  with  Rotational 
•Symmetry  - , - - , - — ~ 


For  a shell  with  rotational  symmetry  about  its  axis,  the 
aerodynamic  force  and  torque  components  will  not  depend  on 
ty" . Change  of  sign  of  6,  A , andZJ,  equivalent  to  increase 
of  0 by  n and  reduction  of  y'by  n,  must  change  the  signs  of 
L,  K,  J,  and  M,  leaving  D and  I unaltered;  while  change  in 

sign  of  £ and  jjF,  equivalent  to  inversion  in  the  origin,  must 
change  the  signs  of  K,  I}  and  J,  leaving  M,  D,  and  L unaltered. 
Thus  if  any  of  these  components  are  split  up  into  terms  even 

and  odd  in  A and  S , the  evenness  or  oddness  of  each  term  in 

and  6 is  determined,  and  we  may  write: 

D » 


VW  pv2  d2  Kj,  ♦ P v d3  sin  6 ♦ p-  ik  sin  6 f £ K„s 

+ . d5  A T,  # v 

’SMS  v B £ A bAS 

L » 

Kl,Kla  pv2  d2  sin  6^  + pv  d3  + P dk  f-  f 


A E Jt 


+ p ~ sin  6 f-  f- 


K = 


kk,kka 


pvd33ln6f.KK  + p<l-f.  TEKKE 

♦ p dU  sin  6 §-  f Kkas 


V ^ „d3 


Via 


i = 

p v d*1  Kj  + p d^  sin  3 + p ^ v sin  6 


+ p d 5 4 f-  KIA£ 


KI2  ,kias 


Vkja 


kj£'kjae 


j * 

p v dU  sin  6 f Kj  + p d5  A £kja  * p dlt  v f KJS 


B A “JA 


+ p d*  sin  6 4-  tt  K 


irr  “JAZ 


M - 


.2  ^3  _4 — v tpvd^f  SiA  + P ^-IZ 


P v d sin  g Kjj  + p 

^•IS^-tAZ 


6a.£1 

v n 13  a 


+ sin6f-f-f-K 


MAZ 


kxf5  ks 

K3Q’,KH 


where  the  coefficient  Kp,  — — ^KA£J  are  •^unctions  M ^d 

R,  even  in  ~ ^ ^ ^ , and  6,  and  are  thereby  uniquely 

determined.  (Note  that  use  of  sin  6 rather  than  6 all )ws  6 to 
pass  through  the  value  n as  well  as  the  value  zero,  our  state- 
ments remaining  validj  but  it  would  not  be  true  to  s y 'hat 
the  coefficients  are  functions  of  sin  6 as  they  will  usually 
differ  at  6 = 0 and  6 » n). 


Tw’'.or1  when  sin  6 ■ 0,  D and  I can  depend  onXand^d  only 

2 2 2 
thr  ough  £ + A , so  that  Kp  and  vail  depend  onS  + 

A c and  *nd  Kj^vnLll  vanish,  while  L and  -K,  and  J and 

M must  behave  like  the  components  of  plane  vectors  depending 
on  the  vector  with  components  £ and  A , so  that  = 

_1Sca  3 hA  = KK£’  ' KJA  * Kime  ’ and  KJS°  *UA>  being  whal 
are  commonly  called  K^j,,  Kg,  K^,  and  K^,  and  likewise  must 


2 2 

be  functions  of57  + A only.  In  general  these  are  all  dif- 

2 2 

ferent  functions  of  £ , + A for  6*0  and  for  6 « n. 


It  must  be  emphasized,  however,  that  when  sin  6 / 0,  these 
last  relations  do  not  havq  to  hold,  although  the  departure  from 
them  must  be  of  order  sin^  6. 
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Part  III.  Motion  of  a Spinning  Symmetrical 
Projectile  with  Large  Yaw 


3.1  The  General  Method  of  Transformation  by  Variation  of  Parameters 

A solution  is  first  made  of  a reduced  problem  in  which  many  terms 
in  the  rates  of  change  of  the  variables  have  been  left  out.  This  leads 
to  a transformation  from  the  variables  to  a new  set,  the  parameters, 
which  would  be  constant  in  time  if  the  retained  terms  alone  existed. 
Applying  this  transformation  to  the  original  equations,  new  exact  equa- 
tions are  obtained  for  the  rates  of  change  of  the  parameters,  which  change 
more  slowly  than  the  original  variables. 

If  both  the  original  equations  and  the  reduced  equations  are  of 
Hamiltonian  form,  the  transformation  may  be  carried  out  on  the  Hamilton- 
ian function  only.  If  the  reduced  equations  are  of  Hamiltonian  form, 
it  may  still  be_  convenient  to  solve  them  by  the  Hamilton-Jacobi  method, 
although  the  transformations  resulting  have  to  be  carried  out  on  the 
separate  equations  instead  of  on  a single  function. 

If  the  reduced  equations  can  be  taken  linear,  we  have  the  usual 
method  of  variation  of  parameters,  which,  for  a second  order  system, 
leads  to  the  W»  K.  B.  method  so  much  used  in  quantum  mechanics. 

If  the  solution  of  the  reduced  problem  is  periodic  with  a compara- 
tively short  period,  the  next  step  usual  in  the  solution  is  to  approxi- 
mate the  equations  of  change  of  the  parameters  by  averaging  over  a period. 
This  gives  an  approximate  system  of  'secular  equations'  for  the  change 
of  the  parameters,  which  represent  the  motion  well,  unless  there  is 
resonance  “with  some  of  the  ignored  periodic  terms. 

If  in  this  last  case  the  reduced  equations  are  in  Hamiltonian 
form,  the  general  increase  or  decrease  in  amplitude  of  the  periodic 
motion  may  be  tested  by  the  average  change  in  the  Hamiltonian  function 
itself,  the  standard  'energy  test  of  stability*. 

It  is  this  last  case  that  we  are  going  to  consider;  so  we  are 
ruling  out  possible  cases  where  very  large  damping  may  make  stable  motion 
that  would  be  unstable  for  small  damping. 

3.2  The  Reduced  Equations  and  Their  Solution 

The  reduced  equations 

6 * A 

0 .*  $ - cos  6 



B sin  6 
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cos  5 
8 sin2  3 


( <f>  - cos  6) 


A = 


- ( § 


cos  6)  ( tz  £ cos  5) 
B sin^  5 


S V(6) 

•56 


^ a 0 

it*  o 

give  the  motion  in  angle  when  all  the  force  and  torque  com- 
ponents except  the  overturning  moment  (K  or  G.)  are  neglected, 
V(6)  and  when  that  is  regarded  as  derived  from  a potential  V(S). 
These  equations  correspond  to  a Hamiltonian  function 


H 


H * 


|B  |-A2  (#  -^cos  6)2j  + |Af 2 + V(6)  , 


just  that  for  a symmetrical  top  with  potential  energy  V(6), 

The  Hamilton- Jacobi  equation 

2 . .2 


as 

"St 


+ i f /as\  + l /as  - as  R\  z\  .1  /as\ 

■sb | (srJ  ^r5\rpf  af  003  6J  J +jk  Isyv 


+ V(6)  * 0 


has  the  complete  integral 
S * S - Ht 

n "* 


t * "ft  f*  f I2®  ( H -JJ-  ~ v(6)) 

° \ 1 i ° 

- (#0-  Jocose)2]7'  Ijf-*, 

£o,fo  and  the  solutions  of  the  reduced  equations  are  given  by 

0 

where  0q,  and  6q,  as  well  as  and  H 0>  correspond  to  time  tQ. 
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5 

5- 


Since  ^ and  tjr  do  not  occur  explicitly  in  H,  £ and  J f are 
equal  to  their  initial  values  $ and  jP  . The  motion  in  6 is 
in  general  a libration  between  two  simple  roots  of 


between  which  this  expression  is  positive. 

It  is  convenient  to  take  6Q  to  be  the  smaller  of  these  roots, 
simplifying  the  partial  derivatives  of  S.  We  call  the  larger  root 


The  solution  is  then  given  by 

1 

A *■  -j  2B  ^ f2  - V(6))  sin2  6 - ( -j^cos  6)2j? 

sin  6. 


§-  § 
tf-  t 


t * t -* 
o 


I 


B sin  6 d 6 


*«■'-/  7 7 — 

6°{2B  (h-I 


2B  (H  - \^fQ  - V(8)  sin2  6 ^cos  6) 

(#0  “ to  003  6) 


(h  - 1st2  - v<6))  sin2  8 - cos  6)’ 


,^1  sin  6 


- # sin41  6 - (#  -3?  cos  6)  cos  6 
A r o '■‘■0*0 


.2 


t-t- 


d5 


2B  (h  -Jif*  “ V(6)^sin26  - ($o~f0  cos  6)2 


1 1 sin  6 


since  the  terms  in  d 6 vanish. 

o 


3.3  The  Uniformising  Variable 

For  manipulating  the  above  integrals  it  is  convenient  to  in- 
troduce a uniformising  variable  u which  continually  increases  as  6 
librates.  Using  co3  6 rather  than  6 as  the  variable  to  be  given  in 
terms  of  u,  we  put 

2 1 2 1 
cos  6 ■ cos  6q  cos  -jju  + cos  6^  sin  -ju 
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u 


sc>  that  6 ■ 6q  for  u * 0,  6 = 6-^  for  u * n,  and  so  on. 

Then  u is  given  in  terms  of  6 by 


tan 


u 


cos  - cos  6 
o 

cos  6 - cos  6, 


K(6) 


with  positive  sign  when  6 is  increasing,  negative  sign  when  6 is 
diminishing. 


Thus 


and  if 


sin  8 d 6 


=*  du 


{” 


■j(cos  6q  - cos  6)  (cos  6 - cos  6ijj'!5 

*» 

(H  • H fo  2 - V(6))  sin2  6 " <*o  -to  003  6)2  ’ 

(cos  6q  - cos  6) (cos  6 - cos  6^)  K(6) 


so  that,  by  the  definitions  of  oQ  and  6^,  K( 6)  does  not  vanish 

for  6 $6^6,,  we  have 

o 1 


% ■ * -Jl  ■yfer 

o 

The  standard  astronomical  proceeding  is  to  expand  1/  / K(6) in  a 
Fourier  series  in  u and  so  obtain  a series  which  can  be  inverted  to 
give  u in  terms  of  t,  although  most  computations  may  be  more  easily 
carried  out  in  terms  of  the  intermediate  variable  u. 

If  we  write  | ^ * 6^  * A + £ , so  that 


tan 


u „/cos 
1?  / cos 


( \ - g ) - cos  6 
6 - cos  (\  + €) 
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we  see  that  u is  the  angle  includeu^W^des  « and  X.of  a spherical 
triangle  with  sides  c , X , and (See' appendix,,) 


Tifi+  tan  v _/cos  (S  - X ) - cos  \ 

‘S  l/  COS  & “ cos  (6  + X ) 


so  that  v and  w are  the  angles  included  by  X and  6,  and  by  6 and  € 
while 

du  sin  6 

w t “S  1 

« (cos  ( 4 -^)  - cos  6)  (cos  6 - cos  ( C +/)) 

dv  _ cos  6 cos  € - cos  X 

cRT*  " r M 

sin  6 •<  (cos  (€  - X)  - cos  6) (cos  8 - cos  (€.+  X))|S 

dw  cos  6 cos  X - cos  £ 

ar"  ? r, 

sin  6 ^ (cos  (€  - X)  - cos  6)(cos  6 - cos  ( * + X)H  ‘ 


so  that 


s X - fo  cos  f)  dv  - (<fo  cos  € - jj^o  cos  X)  dw 


(cos2  € - cos2  X)  / K(6) 

du  + ( cos  X - <ji  Q cos  e)  dv  - <?0  cos  tf-^cosX)  ^w 


The  square  roots  in  the  expressions  for  tan  ^ , tan  -g  , 
w du.  dv  dw 

tan  , and  for  ^ , and  , are  to  be  taken  positive 

when  6 is  increasing,  negative  when  6 is  diminishing,  and  i/"k 
is  to  be  taken  always  positive. 


There  are  now  two  cases 


As  6 increases  from  X - £ to  X + u increases  from  0 
to  Ttj  v increases  from  0 to  a maximum  and  falls  back  to  0,  and 
w decreases  from  n to  0;  as  5 decreases  from  X +€to  X -€  again, 
u increases  from  n to  2 n,  v decreases  from  0 to  a minimum  and 
rises  back  to  0,  and  w decreases  from  0 to  - n;  and  so  on,  u, 
the  uniformising  variable , continually  increases,  v oscillates, 
and  w continually  decreases. 

(b)  <r  > x 


As  6 increases  from  « - X to  X + <,  u increases  from 
0 to  n,  v decreases  from  ji  to  0,  and  w increases  from  0 to  a 
maximum  and  fall3  back  to  Oj  a3  6 decreases  from  X +€to  € - X 
again,  u increases  from  n to  2n,  v decreases  from  0 to  - n, 
and  w decreases  from  0 to  ~n,  and  w decreases  from  0 to  a 
minimum  and  rises  back  to  0;  and  so  on,  u continually  increase, 
v continually  decreases,  and  w oscillates. 


In  the  boundary  case  of  f “ X,  6 passes  through  the  value 
0.  As  6 increases  from  0 to  € + X,  u increases  from  0 to  n, 
v decreases  from  to  0,  and  w decreases  from  ^ to  Oj  as  & 

decreases  from£+  X to  0 again,  u increases  from  n to  2n,  v 
decreases  from  0 to  - ^ , and  w decreases  from  0 to  ; 

as  6 passed  through  the  value  0,  u continues  to  increase,  and 
v and  w may  be  supposed  to  jump  to  the  values  ^ and  ^ again, 

without  change  in  those  functions  of  u,  v,  and  w that  give 
direction  cosines  of  axes  in  the  projectile. 
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Physically,  X is  a mean  yaw,  e an  amplitude  of  nutation. 


3.U  The  Transformation  to  the  Parameters 

VJe  now  use  the  above  equations  to  give  a transformation  from 
the  variables  A,  $,  f}  6,  0,  and  t>  to  the  parameters  H, 

f0  > #0>  tQ  , 0Q  and  fQi  or  to  H,  X,  * , tQ,  0 , and  fQ 

evaluated  at  each  value  of  t,  so  that  V may  be  allowed  to  vary 
with  t as  well  as  with  6. 

__  The  equations  __ 

V(6,t)  H - ~ {a2  + ~~  ( § - f cos  6)2|  + iy  i2  + V (6,t) 

‘-'I  sin  6 J 1 

$o  = $ 

f0  - f 

show  that  we  need  not  distinguish  between  cf  q and  g5  and  between  'Tp” 

andiP,  and,  H - S-  AA  + (#  - J?  cos  6)(  $ - cos  b) 

- x B sin  6 

($-?cos  -e)2»  (f  $ooo  6)?sin  6 6 + §§  s *||  . 


“■f  (°3  ^ I)  ' bW6  < 5 ' 3?  cos  6>  °B  * f Si-ls  K + ft 

F nT  F 

while  $ » GA  cos  6 - G0  sin  6 + A ^ + i — H§ 

and  3^“  = G. 


while 


( § - f cos  5 )(  £ - ■$  cos  6)  + ^ ^ 

B sin  6 

o 


f < Si  (^^C0S  ^ * tSi(^C°3  S“^8ln  6°  f^)a  F 


1*  1 —f 


and  a corresponding  equation  for  6^  give  us  bQ  and  6^  and  so 
X and  k , 


Tne  equation 


where 


t ■ t 
o 


T Bdu 

~7m 


£ 2B  f2  - V (6,  t)  sin2  6 - ( # - $ cos  6)2j-  = 


(cos  6q  - cos  5)(cos  6 - cos  6^)  K(6) 


gives  us 


*0-  1- 


/ K(6) 


u 

“X 


B 


2(K(6)) 


3/2 


6T  tK(6>]  du 


where 


6_ 

6t 


[k(6>]  6F  Lk(6)J  means  lr~  “Hr*  + 


9K  i . 8K  • 
& ds0  6o 


SK  6,  . dK 
3BT  1 + ‘"5h“ 


evaluated  after  6 has  been  expressed  in  terns  of  6q,  6-^,  and  u. 


* B 

Hence  t = 

° / K(6) 

f du 

( w 

fk  3u  j.  _ du 

mv  "SB  °o  BB, 

0 1 

K ) * 

u 

r 

B 6 

[K(6)J  du 

4, 

2(K(6))^  6t 

where  although  the  separate  terns  of  the  first  bracket  become  in- 
finitely large  as  6— ►6Q  or  6— 6p  their  sum  remains  finite 

and  small. 
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• * 

In  a similar  way  0Q  and  t/Tq  can  be  found. 

Thus  the  time  rates  of  change  of  H,  § t f > V 0O,  and  f o, 

3V 

are  found  as  linear  combinations  of  Fg,  F^j  + "gg — > 

3V  t 

and  , with  coefficients  that  involve  § J W j V 0Q,  and 

and  t.  If  these  equations,  and  those  giving  x,  y,  z,  v,  a,  and  0, 
are  averaged  over  a period  of  the  solution  of  the  reduced  equations, 
we  get  the  secular  equations  that  give  the  first  approximation  to-" 

the  changes  of  H,  ^ tQ,  0q,  and  T/Tq.  If  vie  have  been  able 

to  choose  V so  that  tlie3e  changes  are  small  over  a period,  this  will 
be  a useful  method  of  solution.  For  a symmetrical  projectile  the 
right  hand  sides  of  the  secular  equations  will  not  involve  t and  ifT 

and  will  involve  0 only  in  the  terms  involving  gravitational 
forces.  ° 

3.5  Stability  of  Yawing  Motion 


The  kind  of  stability  we  are  concerned  with  is  that  6 and 
6^  should,  if  small,  remain  small,  if  large,  become  small, ’^regard- 
ing^ all  the  variables  except  0,  6,  , and  A,  as  practically 

3TT 

constant.  Taking  V(6)  so  that  + G^  is  small,  we  look  first 

at  reduced  equations  which  give  a conservative  system  with  H and  $ 
constant. 

In  order  that  6q  and  6^  can  be  small,  it  is  necessary  that 

for  the  value  of  $ for  which  £ can  be  small,  namely  the  part 
of  H depending  on  6,  which  reduces  to 


1 — COS  6 1 ir/cN 

3r  rrjsre + v(6) 

should  have  a minimum  at  6 ■ 0.  Then  6^  must  initially  be  less 

than  the  first  maximum  of  this  same  expression,  if  there  is  to  be 
stability. 

• •. 

Next  we  must  look  at  the  secular  equations  for  H and  §>  • These, 
if  the  motion  with  small  yaw  is  stable,  give,  for  some  region  near 
H * V(0),  # - # , motion  tending  asymptotically  to  these  values. 
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This  region  is  the  region  of  initial  values  for  which  the  yawing 
motion  is  to  be  regarded  as  stable. 

It  should  be  remarked  that  this  'secular  stability1  does  not 
exclude  the  possibility  that  in  higher  order  approximation  resonance 
may  exist  leading  to  unstable  motions. 

To  determine  the  region  of  stability  we  solve  the  ordinary 
differential  system  giving  the  secular  change  in  H and  $> , 

H - f(H,  $ ) 

i - g(H,  0 ) 

starting  in  various  directions  from  points  near  the  point  of 

equilibrium  H ■ V(o),  £ ■ , and  going  backwards  in  time,  and  these 

solutions  mark  out  the  region. 

Some  information  may  be  obtained  by  examining  the  equations 

f(H,  £>  ) - o 

g(H,  # ) - o 

for  other  points  of  equilibrium,  but  there  is  no  easier  way  of  deter- 
mining the  existence  of  an  asymptotic  periodic  solution  that  may 
separate  the  stable  and  unstable  regions  than  by  making  a numerical 
integration.  Such  as  integration  for  a single  starting  value  may  be 
sufficient. 

Stable  periodic  solutions  may  also  exist. 

We  have,  in  fact,  to  deal  with  the  simplest  typical  problem  of 
non-linear  dynamics. 

If  we  wish  to  take  into  account  the  effects  of  slow  changes  in  the 

coefficients  and  in  ^ , we  should  work  rather  in  terms  of  X and  € than 

in  terms  of  H and  $ , since  it  is  the  size  of  X and  € in  which  we  are 
really  interested. 
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Part  IV.  Detailed  Approximate  Theory 


U.l  A Special  Form  of  the  Reduced  Equations 

We  may  choose  such  a form  for  V(6)  that  VK(6)  shall  be  constant,  /I 

say. 

Thus  V(6)  - .H  t |B  W„2  -■p/l2 

r 2 c521]r2l/5  2 

M + I jQj  (cos^l  + COS2€  ) - o - o J*  Z\  ° ^ 0 eosXcos  r)cos  6 

2 B sin2  6 


If,  then  we  take 


V(6)  - V 


+ 


2 2 2 
we  have  A (cos  X - cos  € ) “ ( y **  S?  ) +16  Bv, 


p a ~ 

(cos  X + cos  * ) a Q + ) - 16  Bn, 


n 2 n 2 , 2 

a„d“-  -V,  . 


The  solution  of  the  reduced  equations  is  now,  for  € < X, 
^ » «Tl  sin  X sin  v 

# ■ £ o 

V- V. 


t - 1 - Bu/n 


- <p  - i o cos  X - o cos  6 v - (*Q  cos  € - IP  o cos  X)  (n  - w) 

0 (cos^  £ - cos^  X)  (cos^  6 - cos2  X)  ^ 


(* ' x)  re  11  fy°  e°|  a - 

\ cos  € - i 


%cos  i)h  - 

cos c X / 


0 COS  € ~ ^0  COS  X \ ( Tt  - w). 

COS^  € - cos^  X ^ 
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Notice  that  for  V(6)  constant,  this  reduces  to  the  usual  Eulerian 
nutation  of  a symmetrical  rotator,  while  by  taking  v * o,  we  obtain 
motion  corresponding  to  overturning  moment 


k p,  sin  6 
(1  + cos  6)^ 

which  will  serve  as  a good  reduced  motion  for  investigating  stability 
with  large  yaw. 


.In  general,  we  have  a motion  with  mean  yaw  X,  amplitude  of 
nutation  £ and  frequency  of  nutation  «fl/  2 n B. 

The  effects  of  replacing  the  overturning  moment  by  — n — 2 

(l  + cos  6) 

vri.ll  not  show  up  in  our  approximate  system  of#  secular  equations; 

they  show  up  immediately  in  the  equation  for  t which  now  takes  the  form 

*0  ■ TTiS  X Sir:  V ( 4 * i c°3  w + x °°s  v ) > 

containing,  when  € and  X are  expressed  in  terms  of  the  forces,  the 


term 


B 

SV2 


(cot  € cos  w + cot  X cos  v)  . ( 0,  - lfe.5tB.-Q 

(l  + cos  6)‘ 


2 3 

The  leading  term  in  * M is  p v d sin  6 , 

and  we  should  take  ft  so  that 

(cot  € cos  w + cot  X cos  v)  (p  v2d^  sin  6 K„  - - , 

(l  + cos  6)  ' 

vanishes  when  averaged  over  the  reduced  motion,  after  inserting  the 
actual  variation  of  with  6. 

If  does  not  depend  on  6,  this  gives  (see  Appendix) 

2 cos  € cos  X p v2  d^  L = -r 

n cos  c + cos  X 

In  this  case  x,i"5‘5‘l~F*‘  + ^as  no  max^jnuin#  and  ^ae  only 

2 

preliminary  condition  is  IP"  ? 1;  B ft.  that  is,  that  the  stability 
factor  for  small  yaw  is  greater  than  unity.  If,  however, 

p 

increases  with  6 faster  than  1/(1  + cos  6)  , a finite  stable  range 

of  6 may  be  determined. 
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4.2  An  Approximate  Form  for  the  Secular  Equations 


We  shall  systematically  neglect  terms  of  higher  order  than 
the  first  in  A and £in  the  aerodynamic  force  and  torque  components, 
so  that,  for  instance,  we  take 

J - p v d^  sin  6 ^ Kj  - p d^  ^ ^ + p d*  , 


with  the  coefficients  taken  to  be  functions  of  M and  E,  and  everr- 

“<r 

v A 


functions  of  ■-  ? and  8,  but  not  to  depend  on  A and  H . 


We  shall  also  omit  the  gravitational  terms. 

Since  in  the  secular  equations  terms  linear  in  A will  vanish 
on  the  average,  we  obtain  

A V A.  - jb 


H 


,4  /a\2 
•pvd  - 1 

\ B/ 

1 K + o v > V % . 

Kjj  + p V d -J  B ! 

u(ff 

Kj.  -p  vd^ffsin 

-pvd  LA y 
4 T V 

_P_£  ff 

Mv  li  B 

3v  v 3v 

35  KXF  + 3F 

,4 

“-pvd  ■¥ 

cos  S L -pvd 

2 2 
v d 


(S'? 

V / 


% 


dv 
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*L 


P v d*4  % sin2  6 Kj  -pvd1*  Ig  sin  6 - § cos  6^  + 


P d^  IP  (j]f-  <$>  cos  s')  S'  v n V a2 
-JftT  J ST 7TT. U Tf  ^yt?  + P 4 ■ JA. 


sin  6 


S XF 


lv"  ITkxf 


and 

¥- 


-P  vdu  f-ic- 


- p v ^ sin  6 ^ 


These  averages  should  be  taken  over  the  reduced  motion  after 
inserting  the  actual  variations  'of  the  coefficients  with  6.  This 
may  perhaps  be  done  by  expressing  them  approximately  as  polynomials 
in  6. 

• • • 

We  will  now  have  H,  § , and  , given  in  terms  of  , Tp , X,  € , 
and  .ft  , as  well  as  £,  and  we  can  use  the  equations 


cos  X - cos  4 
cos  X + cos  4 


- ( g - ff)  /A 

■/(  $ )2  - 16  b » /n. 
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where  fl  - Ji  B (H  - -p  *tfr2  - VQ) 

cl  d 

to  obtain  ^ cos  X and  ^ cos  t in  terms  of  X,  c , Ip  , and  |4.  The 
stability  problem  then  involves  the  corresponding  problem  in  non-linear 
dynamics,  regarding  p and  ^ as  slowly  changing. 

If  the  aerodynamic  coefficients  do  not  depend  on  6 we  obtain 
(see  Appendix) 


H 


P v d^ifl2  sin2  X 1 - cos  £ 

~ 5— 

B sin*  X 


+ p v d 


^r^  + p cose  cos  X - ^ Kj 

- cos  X) 


h W f cos  € \ 2 cos  X A $ 2 c 

— . — »-= — - — 5 1 “ ~T  5 2 “2  T 

u cos* € - cos c X / B cos  e - cos*  X Bd  cos  £ 


- p v d1 


KI  - p ’ dU  ^?«»«««  * - f )KIZ*  W1-  „ — -1)  *1 


* „ dh  v ,../V  COS2£  + cos2  X + COS£  coo  X-l  $ 1 + cos  £ cos  X\t 

P^jl,V  (=oSf'~07t?  ^ (oos^cosxp)' 


XF 


cos  £ +COSX 


- p v d^  cose  cos  X Kj.  - p v d^^^“(cos2£  cos2  X + -^sin2*  sin2X^  - ^-cos£  cos^K^^j 

- p V d^  Y ^l-cos2£  cos2X  - ■jsin2  £ si.2X  ^ Kj  - p v d^^  cos  £ cos  ^ ^ ^ % 


2 2 

+ -#cos€cos  X^K^  + 


P d^  iP'J?2  +^2  cos  } 

*r"T~ 


o 

cos  X 


$7?  / 2 cos  6 

TT  ( 


\ 2 2 
\ cos  £ -s  cos  X 

+ £^L  f.  fi  six?  \ 

-W-  * —g 


-1) 


1 - COS  £ 

• 2 - 
sxn  X 

$r»-pvdli  f.  Kj  - p v d^  £■ 


K 


XF 


cos  € cos  X - 


# 


) KI 2 
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U.3  The  Special  Case  of  Small  faw 

1 1 ~ ' 1 * 1 O 0 * * 

The  equation  A2  (cos  X - cos  6 ) - ensures  t-hsi  H,  $ , 

and  Tjr.  are  not  large,  on  account  of  the  denominator  cos  X - cos 6 , 
when  the  yaw  is  small.  Using  this  equation  and  keeping  only  terms 

of  order  - § and  £ 2 and  X2  we  obtain 

» = o v d ufa2  \l  *€2  * -ti*-  £i— ]k 

\T  ~ ^ B2  B2  J* 

+ p v d1*  - f f\ 

-pvd“ 

* P 4 ! > (V1  - ? - 4 + ^-) 

, du  | £ . - p vd“  [stf*l  - f (^|IS 


.„d^|  u’u’iyp^  < % 

. p ^ V-  t + £sf-  ¥•  iar  (‘2+  x2)  'f  (_t:LX)kj 

^..p.d1*  f Kj  -pvd^lii  -f  i-(£2»X2)]KI2 


For  small  yaw,  X and  € are  small,  and 

A2^T  2 -UBti 


? 3 

p,  *=»  p v dJ  Kjj 

^A2  .2 


while 


Tfr  - $ 7^^^- 


■gr2  0 „ A2  m2..  2, 

and  H - + 2p,  - VQ  ^ -jjg-  (x  + 6 / 


I 


3U 


n 


o 2 

These  are  of  the  first  order  in  X and.  £ , and  so  are  their  time 

rates  of  change. 

In  fact,  if  we  write  for  brevity, 


h-“£JL£ 

B 


p v d 


Mv 


L P v d 

* ~T“ 


*l 


Y « P.JLd. 
> A 

✓v  P d^ 


XF 


we  obtain 

¥ 


jr_^=^|^Y+^h-Jk-|^  + n't  J ^X2  + £ 2)-J^  h + ^ k)(x2"  € 2 ) 

h - +2 i - ^ k(ip2  -A2)-  ^ h($2  + A2)  -tf2( Y + H*) 

- B njj-^X2  + * 2)  + (h  * Y + ‘X  (x2  - € 2^  , 


while,  keeping  only^the  leading  term, 

' . A 

A cr  d L T7)r 

From  these  equations  we  find,  writing  a * , — * j>y~  “ ~^r 


so  that 


-f -J(,  "hi1®  Bi-  ^ A2  (i-|)  . 

k if  - * 2 ) • 5 C2r  ♦ » - 11  - ■*  * 2 * I*)^2  * <2)  '(h  * k * s)p-  ‘ 2) 

4 (X2  . * 2)  - -('h  ♦ k -i  . ?J(x2  .i  ^2T  + h - k - 2Xk)^2  -f2) 


so  that 


d ,2 


cTE 


Xe  » 


- ♦ k - j?  + | ) - i (by  ♦ h - k -jf  ♦ 2*  n) 


J 
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and  the  condition  of  stability,  for  real  a,  is  that  both  the  coefficients 
on  the  right  hand  side  should  be  positive. 


APPENDIX 


Formulas  connected  with  a spherical  triangle 


X 

If  u5  v,  and  w,  are  the  angles  of  a spherical  triangle  with 
sides  6,  € , and  X,  we  have  the  following  formulas. 

sin  u a sin  v a sin  w 
sin  6 sin  6 sin  X 


tan 

u 

. /cos 

(X  - c ) “ cos 

1 

6 

2 ' 

V cos 

5 - cos  (X  + € 

7 

tan 

v /cos 

(6  ~ X)  ~ cos  £ 

Y cos 

€ - cos  (6  + X) 

tan 

w / cos 

(£  - 6)  ~ cos  X 

V cos 

X - cos  (.£  + 6) 

cos 

6 

= cos  6 

cos  X + sin  fe  sin  X 

cos  u 

cos 

e 

* cos  X 

cos  6 + sin  X sin  6 

COS  V 

cos 

X 

* cos  6 

cos  « 

+ sin  6 sin  £ 

cos  w 

cos 

u 

m ~ 

cos 

v cos  w + sin  v 3in  w 

cos  6 

cos 

V 

m — 

cos 

w cos  u + sin  w sin  u 

cos  € 

cos 

w 

x — 

cos 

u cos  v + sin  u sin  v 

cos  X 

cos 

V 

si;: 

6 = 

sin  X 

cos  £ - cos  X 

sin  £ cos  u 

cos 

V 

sin 

X * 

sin  6 

C03C  - cos  6 

sin  € cos  w 

cos 

w 

sin 

6 - 

3in  £ 

COS  X - COS  C 

sin  X cos  u 

cos 

w 

sin 

€ * 

sin  6 

cos  X - cos  6 

sin  X 

cos  V 

cos 

u 

sin 

X * 

sin  € 

COS  6 - C06  £ 

sin  6 

cos  w 

cos 

u 

sin 

€ « 

sin  X 

cos  6 - cos  X 

sin  6 

cos  V 

cos 

fe 

sin 

w * 

sin  u 

cos  V + cos  u 

sin  v 

cos  X 

cos 

6 

sin 

u * 

sin  w 

COS  V + cos  w 

sin  v 

cos  6 
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cos  X 

sin 

u 

CB 

sin  v 

cos  w + 

cos 

V 

3in  w 

cos  6 

cos 

X 

sin 

V 

IS 

sin  u 

COO  W + 

cos 

u 

sin  w 

cos  «. 

cos 

S 

sin 

V 

M 

sin  w 

cos  u + 

cos 

w 

sin  u 

cos  6 

cos 

6 

sin  w 

M 

sin  v 

cos  u + 

cos 

V 

sin  u 

cos  X 

• 

• 

• 

• 

u 

sin 

V 

sin  X * 8 

- 6 cos 

w - 

X 

cos  V 

* 

• 

* 

•» 

V 

sin 

w 

sin  8 = £ 

- X cos  u - 6 

cos  w 

* 

w 

sin 

u 

• 

sin  £ ■ X 

• 

- 6 cos 

v «» 

# 

£ 

cos  u 

The  above  formulas  mav  be  used  to  average  expressions  with 
respect  to  u.  In  particular,  since 


du 

a + b cos  u 


* 


and 


we  find,  for 


t 2n 
*/ 


du 

— ? 

(a  + b cos  u)‘ 


£<X,  the  following  averages: 


3/2 

) 


cos 


* cos  € cosX 


p 2 12  2 

2C  * cos  € cos  X + sin  € sin  X 

COS  O £ 


1 

1 + cos  5 


COS  € + COS  X 


1 - cos  6 


COS  € - COS  X 


sin^  6 


■ C03-f 

.2  2 


COS  e.  - cos  X 


cos  6 


‘A 

sin*  6 


cos  € - cos  X 
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